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' The well known general relativistic Lense-Thirring drag of the orbit of a test particle 

in the stationary field of a central slowly rotating body is generated, in the weak-field 
and slow-motion approximation of General Relativity, by a gravitomagnetic Lorentz-like 
acceleration in the equations of motion of the test particle. In it the gravitomagnetic 
field is due to the central body's angular momentum supposed to be constant. In the 
context of the gravitational analogue of the Larmor theorem, such acceleration looks like 
a Coriolis inertial term in an accelerated frame. In this paper the effect of the variation in 
time of the central body's angular momentum on the orbit of a test mass is considered. 
O . it can be shown that it is analogue to the inertial acceleration due to the time derivative 

of the angular velocity vector of an accelerated frame. The possibility of detecting such 
> * j ' effect in the gravitational field of the Earth with LAGEOS-like satellites is investigated. 

j^JT), It turns out that the orbital effects are far too small to be measured. 
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1. Introduction 



In the weak-field and slow-motion approximation of General Relativity, a test par- 
ticle orbiting a central slowly rotating body of mass M and angular momentum J, 
assumed to be constant, is acted upon by a non-central acceleration 1 of the form 

a G M = - x B g , (1) 

in which v is the velocity of the test particle, c is the speed of light in vacuum and 
B g is the gravitomagnetic field given by 2 

2g[J-3(J-f)f] 

c r A 

In it f is the unit position vector of the test particle and G is the Newtonian grav- 
itational constant. The gravitomagnetic field B g can be derived 2 from the gravito- 
magnetic potential A g 

A g = ^ (3) 



according to 



VxA s = -| (4) 



i 
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Equation (j|) and Equation (^) hold far from the central body supposed to be 
spherically symmetric and rigidly rotating. An inertial frame K{x, y, z} with its 
origin located at the center of mass of the rotating body, the z axis directed along 
J and the {x, y} plane equal to the equatorial plane of the gravitating source is 
adopted. 

This effect, which affects the longitude of the ascending node fl and the argu- 
ment of pericenter oj of the orbit of a test particle 3 , is currently under measurement 
by analyzing the laser-ranged data to LAGEOS satellites in the gravitational field 
of the Earth 4 . 

According to the gravitational analogue of the Larmor theorem 5 , we could ob- 
tain Equation (|l|) by considering an accelerated frame rotating with angular velocity 

= ff. (5) 

Indeed, in it an inertial Coriolis acceleration 

a Cor = 2v x n L (6) 

is experienced by the proof mass^J 

It seems natural to pose the following question. In an accelerated frame, apart 
from the Coriolis and centrifugal inertial forces^], a particle feels an acceleration 

as well if the angular velocity vector J7l is time-dependent. Does a gravitomagnetic 
analogue of such term exist? If so, it should be induced by temporal variations of 
the angular momentum of the central body J. In the case of the Earth we know 
that J ffi changes in time due to the luni-solar torques which generate, among other 
things, the precession of equinoxes 6 . 

In this paper we will try to investigate this feature. The plan of the work is as 
follows. In Section 2 we derive the full expression of the gravitomagnetic acceleration 
experienced by the test mass when the gravitomagnetic potential due to the central 
rotating mass is time-dependent. In Section 3 the orbital effects on the semimajor 
axis, the inclination and the node of the orbit of an Earth satellite are derived and 
the possibility of measuring such effect in the terrestrial field with LAGEOS-likc 
satellites are discussed. Section 4 is devoted to the conclusions. 

In Table 1 the numerical values of the parameters used in the text are quoted. 
For the Earth, the Moon and the Sun reference 7 has been used. 



a It should be noticed that the equivalence is useful only in spatial regions in which the field 
considered is uniform, so that 5 — 4A g = B g X r. 

b The centrifugal force is proportional to the square of S^l, so that the impact of a possible 
gravitomagnetic analogue, proportional to the square of B g , can be neglected because we are in 
the weak-field approximation. 
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Table 1. Parameters used in the text. 



Parameter 


Description 


Value 


Units 


M 


Mass of the Sun 


1.99 x W' J ' J 


g 


do 


Earth-Sun mean distance 


1.496 x 10 13 


cm 


M m 


Mass of the Moon 


7.3483 x 10 25 


g 


dm 


Earth-Moon mean distance 


3.844 x 10 10 


cm 


(C-A)® 


Earth's principal moments of inertia 


2.6249 x 10 42 


g cm 2 


G 


Newtonian constant of gravitation 


6.673 x 1CT 8 


g cm° s 


c 


Speed of light in vacuum 


2.9979 x 10 10 


cm s — 1 


£ 


Obliquity of ecliptic 


23.44 


deg 


fflLAGEOS 


LAGEOS semimajor axis 


1.2270 x 10 9 


cm 


e LAGEOS 


LAGEOS eccentricity 


0.0045 




*LAGEOS 


LAGEOS inclination 


110 


deg 


"LAGEOS 


LAGEOS mean motion 


4.643 x 10" 4 


rad s — 1 



2. The gravitomagnetic acceleration 

In the weak-field and slow-motion linearized approximation of General Relativity, 
the Lagrangian £ of a test particle of mass m moving in a metric ds 2 — goo(dx ) 2 + 
gu{dx 1 ) 2 + 2g kdx°dx k is given by 2 

C = £ (GE) - ^A g • v, (8) 

where the gravitomagnetic potential A g is generated by the off-diagonal components 
gok, k = 1,2, 3 of the space-time metric while £(ge) denotes the gravitoelectric 
Schwarzschild terms. The Lagrange equations of motion are 

d fdC\ DC „ 

From Equation (||) it follows 

d / dC\ dv 2m. . 2mdA„ 

at U = m Tt ~ - (v ' v) A ^ 7¥ + (GE) ' (10) 



dC 2m 
dr c 

Then, the equations of motion of a test particle are 



_ = -_V(A g .v) + (GE). (11) 



dv 2m—,, . 2m, _. . 2mdA„ 
m- = --V (A g • v) + - (v • V) A g + --^ + (GE) 



[vx(VxA s )] + -^|(GE). (12) 



2m , ,_ . 2m 9A„ 

-[vx(VxA g)] + -— ^ 

In Equation (O) (GE) represents the post-Newtonian gravitoelectric acceleration 
due to the Schwarzschild part of the metric, of order 0(c~ 2 ) +O(c ), which reduces 
to the Newtonian central term, of order O(c ), in the limit c — * ao. For — 4A g = 
B g x r, i.e. when the field is uniform, the gravitomagnetic part of Equation ( p^ ) 
becomes 

a G M = -xB g --^xr. (13) 



c 
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If we use Equation ^ , Equation ([T^) becomes 

a GM = 2v x L + r x — (14) 

Equation ( |i"4| ) shows that the gravitational analogue of the Larmor theorem ex- 
tends exactly also to the case in which the gravitomagnetic field is explicitly time- 
dependent. 



3. The consequences of the variability of Earth's angular 
momentum 



We will focus on the second term of Equation (|13|). When the mass-energy distri- 
bution is a central spherically symmetric rigidly rotating body, so that Equation 
(||) and Equation (^) can be applied, Equation ( p"3| ) can be written 

2G dJ 2G 

c 2 r 3 dt 



a L- = -2-3— x r = M x r, (15) 



in which M = ^ is the momentum of the external forces. Indeed, in a body-fixed 
rotating frame with an angular velocity luq the dynamical Eulcr equations hold 

- = -+. xJ = M. (16) 

However, we are using an inertial frame so that 

dJ dJ m _ .„„, 
— = — =M. 17 

dt dt y ' 

In the case of the Earth M is the external torque exerted on its equatorial centrifugal 
bulge by the other bodies of the Solar System 8 . We will consider only the Moon 
and the Sun^. In the frame K previously defined and by assuming the Earth as 
an oblate spheroid whose rotation is affected by the presence of a disturbing body 
of mass Mb (e.g. the Sun) moving around it on an approximately circular orbit of 
radius d lying on the ecliptic plane, the external torque can be written as 8 

M x = 3GMb( ^ ~ A)@ sin 2 A sin e cos e (18) 

3GM B (C-A) B 
My = smAcosAsme (19) 

M z = 0, (20) 

where (C — A)® — J2M (S R^ in which Ji is the first even zonal coefficient of the 
multipolar expansion of the terrestrial gravitational field and M® and Rq are the 

c Recall that the other planets induce secular effects on the ecliptical plane: the planetary precession 
and a secular change in the obliquity e. The luni-solar torque affects the equatorial plane causing 
the Earth's spin to precess and nutate about the ecliptic pole. While the precession does not affect 
e but only the position of the equinox in the ecliptic, the nutation does affect e periodically. 
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mass and the equatorial radius, respectively, of the Earth. Moreover, A is the eclip- 
tical longitude of the perturbing body and e is the inclination of the ecliptic to the 
equator. 

The next steps in order to evaluate the observable consequences on the orbit of 
an Earth artificial satellite consist of projecting the perturbing acceleration onto 
the usual radial, along and cross-track directions, evaluating the so obtained com- 
ponents R, T and N of aLar in terms of the osculating Keplerian elements of the 
orbiter and calculating their rates of changes by averaging over an orbital revolu- 
tion of the satellite. In performing the orbital averages d and e will be considered 
constant. 

From Equation (|l5|) and Equation ([l8])-Equation ( po|) it is straightforward to 
obtain 

R = (21) 
G\ 2 M B {C ~ A) (B smi 



where 



c 



r = e - DV — r (22) 



d 3 r 



X {> (-) 2 MB % r 2 Ah i Acos ^ + /) + 6sin (^ + /)] > ( 23 ) 



c 



_ sin2£sinfi(l — cos2A) sin e cos O sin 2 A , . 

r = i '- + (24) 

4 2 V ' 

sin 2s cos i sin fi (cos 2 A — 1) sin e cos % cos fi sin 2 A . „ 

A= ^ (25) 

„ sin 2e cos O (cos 2 A — 1) cos e sin fl sin 2A , . 

B= 1 " + 2 ' ( 6) 

in which i is the inclination of the satellite's orbit and / is its true anomaly. From 
Equation (pl|) -Equation ( p3| ) it can be noticed that the disturbing acceleration is 
non-central. 

Among the Keplerian orbital elements of a typical Earth satellites the scmimajor 
axis a, the the inclination i and the node f2 are, in general, well measured, as in the 
case of LAGEOS laser-ranged satellites. Then, we will calculate the rates of change 
of such elements according to 6 



da 2 



df nVl 



Re sin(w + /) + T 



a(l - e 2 ) 



(27) 



%= i 2 iV-cos(c + /), (28) 
at naVl - e 2 a 

c\Q A r 

~~tt = 7= 2 . N- sin( W + /), (29) 
oz nay/ 1 — e A smi a 

in which e is the satellite's orbital eccentricity, n = %/ GMa~ 3 is the satellite mean 
motion and 

(») 

1 + e cos / 
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on the unperturbed Keplerian ellipse. 

By inserting Equation (|2l])-Equation ( p3| ) in Equation (^(-Equation (^), av- 
eraging over an orbital revolution of the satellite and neglecting terms of order 
0(e), it is possible to obtain the long-term evolutions of a, i and fl 

2 



da 
dt 



M B (C- A)z 



sin i sin e 



1 



cose 



na 2 d 3 
sin(f2 - 2A) 



[cose sin f2- 



1 — cos e 



di_ 

dt 



2 

GA 2 Mb(C — A)® cosisine 



sin(fi + 2A) 



1 



COS £ 



~dl 



^1=6 - 



2 

G x 2 
c 



4na 3 d 3 
sin(fi - 2A) + 
M B {C-A) ( 



: sin fi- 



i sine 



1 



cose 



Ana 3 d 3 sini 
cos{tt - 2A) - 



1 — cos e 
2 

[cos e cos f2 
1 + cos e 



sin(fi + 2A) 



cos(D, + 2A) 



(31) 



(32) 



(33) 



From an inspection of Equation (|3l|)-Equation ( 33 ) it can be noticed that there 
are no secular, linear trends but only long-period harmonic perturbations whose 
frequencies are linear combinations of the longitude of the satellite's ascending 
node f2 and the ecliptical longitude A of the disturbing body B. E.g., for the Sun 
and LAGEOS we have 

P(ft) = 1043.67 days, (34) 
P{fl - 2A) = -221.34 days, (35) 
P(Q + 2A) = 155.42 days. (36) 

Notice also the d~ 3 dependence, typical of a tidal effect. 

Concerning the possibility of measuring such effects in the terrestrial field, unfor- 
tunately there is no hope of detecting them: indeed, they are c -2 effects multiplied 
by G 2 . For LAGEOS the experienced acceleration due to the solar torque is of the 
order of 



61° 



M Q (C-A) e 



1CT 17 cm s" 2 , (37) 

while the present level of sensitivity is 9 10~ 10 cm s~ 2 . Consequently, the ampli- 
tudes of the periodic signals of the semimajor axis, the inclination and the node 
amounts to 10~ 7 cm and 10~ 8 mas, respectively^. If the action of the lunar torque 
is considered, it turns out that 



"Lar 



"Lar 



2.1. 



(38) 



d In integrating Equation (jjl])-Equation ( jj^ ) the corrections AA and Ae due to precession and 
nutation have been neglected. 
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A posteriori, this result justifies the previous choice of neglecting terms of order 
0(e) in the satellite's eccentricity and the contribution of the torques of the other 
planets of the Solar System. 

4. Conclusions 

Motivated by the gravitational analogue of the Larmor theorem, in this paper we 
have investigated some consequences of the variability of the proper angular mo- 
mentum of a central slowly rigidly rotating body due to external torques on the 
gravitomagnetic equations of motion of a test mass in the weak-field and slow- 
motion approximation of General Relativity. The test particle turns out to be af- 
fected by a post-Newtonian 0(c~ 2 ) non-central acceleration analogue to the inertial 
acceleration due to the temporal derivative of the angular velocity vector arising in 
a non-inertial rotating frame. In view of a possible measurement of such effect with 
LAGEOS-like satellites in the gravitational field of the Earth, we have calculated 
the perturbations induced on the semimajor axis, the inclination and the node of an 
orbiting proof mass. We have found no secular terms but only long-term periodic 
effects: their magnitude is far too small to be detected. 
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